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[l] We investigate a coronal mass ejection (CME) propagating toward Earth on 29 March 2011. This 
event is specifically chosen for its predominately northward directed magnetic field, so that the 
influence from the momentum flux onto Earth can be isolated. We focus our study on understanding 
how a small Earth-directed segment propagates. Mass images are created from the white-light cameras 
onboard STEREO which are also converted into mass height-time maps (mass J-maps). The mass 
tracks on these J-maps correspond to the sheath region between the CME and its associated shock 
front as detected by in situ measurements at LI. A time series of mass measurements from the STEREO 
COR-2A instrument is made along the Earth propagation direction. Qualitatively, this mass time series 
shows a remarkable resemblance to the LI in situ density series. The in situ measurements are used as 
inputs into a three-dimensional (3-D) magnetospheric space weather simulation from the Community 
Coordinated Modeling Center. These simulations display a sudden compression of the magnetosphere 
from the large momentum flux at the leading edge of the CME, and predictions are made for the time 
derivative of the magnetic field (dBldt) on the ground. The predicted dB/dt values were then compared 
with the observations from specific equatorially located ground stations and showed notable similarity. 
This study of the momentum of a CME from the Sun down to its influence on magnetic ground stations 
on Earth is presented as a preliminary proof of concept, such that future attempts may try to use 
remote sensing to create density and velocity time series as inputs to magnetospheric simulations. 

Citation: Savani, N. R, A. Vourlidas, A. Pulkkinen, T. Nieves-Chinchilla, B. Lavraud, and M. J. Owens (2013), 
Tracking the momentum flux of a CME and quantifying its influence on geomagnetically induced currents at 
Earth, Space Weather, 11, 245-261, doi:10.1002/swe.20038. 

1. Introduction 

[ 2 ] The influence of a southward orientated mag- 
netic field of the solar wind impinging onto the Earth's 


'University Corporation for Atmospheric Research (UCAR), 
Boulder, Colorado, USA. 

2 NASA Goddard Space Flight Center, Greenbelt, Maryland, USA. 
3 Space Science Division, Naval Research Laboratory, 
Washington, DC, USA. 

4 IACS, The Catholic University of America, Washington, DC, USA. 
5 Institut de Recherche en Astrophysique et Planetologie, 
Universite de Toulouse (UPS), France. 

6 UMR 5277, Centre National de la Recherche Scientifique, Toulouse, 
France. 

7 Space Environment Physics Group, University of Reading, 
Reading, UK. 

Corresponding author: N. P. Savani, Solar Physics Branch, Naval 
Research Laboratory, 4555 Overlook Ave. SW, Washington, DC 20375, 
USA. (neel.savani02@imperial.ac.uk) 

©2013. American Geophysical Union. All Rights Reserved. 
1542-7390/13/10.1002/swe.20038 


magnetosphere is known to be the main driver for 
coupling the solar wind to the terrestrial system [e.g., 
Gonzalez and Mozer, 1974; Bargatze et al., 1985; Tsurutani 
et ah, 1992a]. In this case, the oppositely oriented magnetic 
fields between the solar wind and the magnetosphere are 
more easily able to reconnect and thereby transfer energy, 
mass, and momentum to the Earth. For this reason, sig- 
nificant emphasis on understanding and predicting the 
magnetic field orientation has been pursued to improve 
our space weather predictive capabilities. 

[ 3 ] The role of the solar wind density, N sw , in space 
weather studies is more complex. Statistical studies 
aimed at coupling the solar wind density to geomagnetic 
indices such as Dst show a weak correlation [O'Brien and 
McPherron, 2000]. Others suggest that N sw should be sig- 
nificant in mediating the energy transferred to the mag- 
netosphere [Borovsky et ah, 1998; Thomsen et ah, 1998]. 
The density has been demonstrated to affect the mag- 
netospheric response by way of the solar wind dynamic 
pressure (momentum flux) [e.g., Xie et ah, 2008]. A sharp 
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increase in solar wind dynamic pressure can rapidly com- 
press the Earth's magnetopause causing a large time rate 
of change in the local ground-based magnetic field inten- 
sity. While this effect is not often considered as a measure 
of "geo-effective" events, the space weather consequences 
can be severe. An example of this can be found in the 4 
August 1972 event, where the peak Dst index was mod- 
erate compared to the resulting outage of the AT&T 
telecommunications cable and the geosynchronous satel- 
lite solar cell damage [Anderson et al., 1974; Lanzerotti, 1992; 
Tsurutani et al., 1992b]. 

[ 4 ] When investigating the drivers of space weather at 
Earth and in particular geomagnetically induced currents 
(GICs), it is important to distinguish between the different 
solar wind structures that initiate the chain of events that 
may lead to significant socioeconomic losses. Different 
interplanetary structures have been previously reported to 
produce significantly varying responses in the Earth mag- 
netosphere [e.g., Huttunen et ah, 2008; Miyoshi and Kataoka, 
2005; Borovsky and Denton, 2006; Denton et ah, 2006; Lavraud 
and Borovsky, 2008], and it is therefore reasonable to expect 
different responses for GICs. The two main drivers of 
GICs are from intense magnetic storms associated with 
interplanetary coronal mass ejections (CMEs) and co- 
rotating interaction regions (CIRs). Kataoka and Pulkkinen 
[2008] and Borovsky and Denton [2006] both found that the 
biggest problems for ground-based conducting systems 
were from CME drivers, whereas CIR-driven storms had 
a relatively minor effect. 

[ 5 ] Fast-moving interplanetary CMEs can be considered 
to be made of two primary constituents: (i) the sheath 
region between the shock front and the leading edge of 
the CME and (ii) the ejecta itself. The ejecta is often mag- 
netically dominated (i.e., plasma ft < 1) and displays the 
magnetic properties of a flux rope (FR) when seen in 
situ [Burlaga, 1988]. The sheath region however contains 
denser compressed plasma with a higher value of plasma 
ft. The direction of the magnetic fields within this region 
is variable, but is often found to vary within a 2-D plane 
[Nakagawa et al, 1989; Jones et al, 2002; Kataoka et al, 2005; 
Savani et al, 2011]. The momentum flux, P dyn = N SW V 2 
within the sheath is typically high. Huttunen et al. [2008] 
concluded that as a geo-effective CME propagates over 
the Earth (usually for a period of —24 hours), the most 
intense GIC activity recorded on the ground is most likely 
to occur at the beginning, during the passage of the sheath 
region, and between the shock and the leading edge of the 
magnetic obstacle. 

[6] A high solar wind density can cause a change in the 
compression ratio of the bow shock, as is frequent within 
a geomagnetic storm [Lopez et al, 2004]. For southward 
magnetic field, this serves to increase the sensitivity of 
the predicted effects to N sw [see also Lavraud and Borovsky, 
2008]. Recent studies have considered the effect of pre- 
conditioning the magnetosphere to increase the space 
weather effects. For geomagnetic storms driven by both 
CMEs and CIRs, Lavraud et al [2006] showed that the 
Dst were underpredicted by a model for intervals that 


were preceded by an extended period of northward IMF. 
This is consistent with the hypothesis that a precondi- 
tioned high density plasma sheet will cause a larger than 
expected geomagnetic storm, as also expected from sim- 
ulations [Lavraud and Jordanova, 2007]. The density of the 
Earth's plasma sheet is regulated by the solar wind den- 
sity and lags behind by a few hours [Borovsky et al, 1998]. 
For CME-driven storms, the plasma sheet is generally 
more dense and persists longer than those driven by CIRs 
[Denton et al, 2006]. A magnetohydrodynamic (MHD) sim- 
ulation studied by Siscoe et al. [2002] found that the polar 
cap potential saturation follows a power law with P dyn . This 
then led Xie et al. [2008] to develop a model where the 
magnetosphere can be preconditioned by P dyn . 

[ 7 ] The Dst index is often used as a measure of 
the severity of a geomagnetic storm. This index is an 
average estimate of the global response as determined by 
four equatorially located ground stations. The Dst mea- 
sures the horizontal component of Earth's ground mag- 
netic field and is predominately monitoring the strength 
of the ring current in the magnetosphere [Fukushima and 
Kamide, 1973; Liemohn et al., 2001]. Since a denser and 
cooler plasma sheet produces a stronger ring current, 
CME-driven storms often produce more intense Dst val- 
ues. Also, low latitude (< 56°) aurora, as measured in 
Japan was found to be predominately CME associated 
[Shiokawa et al, 2005]. During a CME -triggered geomag- 
netic storm, the Dst profile often begins with a short dura- 
tion rise which is triggered by the arrival of the CME shock 
front (sudden storm commencement, SSC). The dura- 
tion of the rise phase between the SSC and storm onset 
(SO) often corresponds to the sheath region upstream of 
the magnetic ejecta of the CME. The main phase of the 
geomagnetic storm, which corresponds to the magnetic 
ejecta traversing over the magnetosphere, is where a large 
decrease in the Dst value occurs. However, the phases of 
a geomagnetic storm can be time shifted under conditions 
where a strong southward Bz persists within the sheath 
region. The amplitude of the Dst decrease during the main 
phase is closely related to the southward magnetic field in 
the IMF and its subsequent magnetic reconnection with 
the magnetosphere. 

[8] The momentum flux of a CME, and in particular 
the sheath region, is higher than that of the ambient 
solar wind. This introduces a changing pressure onto the 
magnetosphere, and in the case of a CME-driven shock 
front, the change is abrupt. The subsequent compression 
of the magnetosphere is expected to generate currents 
within the terrestrial system which are then detected on 
the ground as magnetic field fluctuations. As displayed by 
a Dst time profile during a geomagnetic storm, the geo- 
magnetic field undergoes fast changes over a short time 
period. The largest time derivative of the ground magnetic 
field, dBldt, occurs during the initial phase (between SSC 
and SO) or during substorms. As indicated by Faraday's 
law of induction, dBldt is key in estimating the geomag- 
netically induced currents in technological conductor sys- 
tems like power grids [Pulkkinen et al, 2006; Committee On 
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The Societal and Economic Impacts Of Severe Space Weather 
Events, 2008]. Therefore, dBldt is often considered as a rea- 
sonable proxy for localized GIC activity [Viljanen et al, 
2001]. In order to calculate the geo-electric field, additional 
information on the ground conductivity and specifics on 
the technological systems involved are required [Pirjola, 
2002; Pulkkinen et al, 2007]. In this paper, we will focus on 
tracking the effects of CME momentum flux to the proxy 
of GICs, namely dBldt. 

[ 9 ] Significant geomagnetically induced currents 
have been observed to affect technological systems in 
middle- and low-latitude regions [e.g., Ngwira et al, 2008; 
Liu et al, 2009; Watari et al, 2009], and therefore, it is 
becoming increasingly clear that these geomagnetically 
storm-driven GICs are not just a high-latitude phe- 
nomenon [Pulkkinen et al, 2010a]. The extent to which 
GICs can be expressed by dBldt is also not as straight- 
forward as first expected due to the complications of 
estimating the ground conductivity and how this varies 
at different depths. Pulkkinen et al. [2010a] showed that in 
Japan where the tectonic plates are geologically active, the 
local subduction zone is able to affect GICs in an unusual 
way by following the amplitude of the local geomagnetic 
field rather than the time derivative; this highlights the 
complexity of predicting GICs. 

[ 10 ] In this paper, we focus on better character- 
izing the N sw from remote observations as a solar 
wind driver for dBldt of the geomagnetic field. We 
take steps toward creating a time series of the solar 
wind density and show that the momentum flux of 
a CME is an important parameter when predicting 
space weather incidents. We show this by compar- 
ing the ground station data with predicted estimates 
from the space weather modeling framework (SWMF) 
[Tdth et al., 2005]. 

2. Remote-Sensing Observations 

[ 11 ] The STEREO mission, launched in 2006 [Kaiser et al, 
2008], consists of two spacecraft that follow a trajectory 
similar to that of the Earth. As they separate from each 
other at a rate of ss 45° per year, one spacecraft trav- 
els ahead of the Earth (ST-A) while the other lags behind 
(ST-B). Each spacecraft carries the Sun Earth Connec- 
tion Coronal and Heliospheric Investigation [Howard et al, 
2008] imaging package, which contains an Extreme Ultra- 
violet Imager, two coronagraphs (COR-1 and COR-2), and 
the Heliospheric Imager (HI). The HI instrument on each 
STEREO spacecraft is made up of two wide-held visible- 
light imagers, HI-1 and HI-2 [Eyles et al, 2009], The fields 
of view of HI-1 and HI-2 are of 20° and 70° angular extent, 
respectively, and under ordinary operation are nominally 
centered at 13.7° and 53.4° elongation in the ecliptic plane. 
Thus, the ecliptic plane corresponds to a horizontal line 
that runs through the center of the fields of view. Figure 1 
displays the location of the STEREO spacecraft in rela- 
tion to the Sun and Earth on 27 March 2011. The shaded 
regions indicate the held of view for the HI-1 cameras on 



Figure 1. Location of spacecraft during CME propaga- 
tion on the 27 March 2011. The shaded regions locate 
the held of view for the HI-1 cameras on both the ST-A 
and ST-B spacecraft. 

both spacecraft. The direction of propagation for the CME 
nose is estimated from Jmap techniques [e.g., Sheeley et al, 
1999; Rouillard et al, 2008; Davies et al, 2009; Savani et al, 
2009, 2012a] and detailed below in section 2.3. 

2.1. CME Mass Calculation 

[ 12 ] In order to convert the white-light images that con- 
tain photometric information per pixel into mass per pixel, 
we exploit the Thomson scattering properties of electrons 
from the corona and inner heliosphere. The total light 
observed by the coronagraphic and heliospheric imagers 
onboard STEREO is from photospheric photons scattered 
by all the electrons along the entire line of sight (LOS). 
The difference in contributions of each electron along the 
LOS depends on the distance and the Thomson scatter- 
ing mechanism [Billings, 1966; Vourlidas and Howard, 2006]; 
therefore, an estimate of the total number can be gener- 
ated from white-light images [Hayes et al, 2001]. In order 
to estimate the mass of a CME, the brightness contribu- 
tion of the transient must be isolated from the background 
coronal signal. Previous studies have achieved this by sub- 
tracting a suitably chosen pre-event image from an image 
containing the CME [Stewart et al, 1974; Howard et al, 1985; 
Poland et al, 1981]. A comprehensive explanation of the 
procedures required to make estimates of a CME mass 
from the raw telemetry data received on the ground is 
given by Vourlidas et al. [2010]. However for completeness, 
the main key points in the procedure are given below: 

[ 13 ] 1. First, the relevant time stamps for images con- 
taining the CME are chosen as well as a pre-event image. 
The pre-event image is ideally the last possible image 
prior to the CME entering the field of view in order 
to make the appropriate corrections for background. It 
is important to minimize the effects from evolutionary 
changes and solar rotation, and not to include another 
CME or other transient effects in the pre-event. 
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Figure 2. A mass image from the COR-2, HI-1, and HI-2 cameras from ST-A. The pixel 
intensity displays the locations of the highest line of sight mass measurements. The cone 
overplotted on the figure displays the relevant part of the CME that propagates over 
the Earth. 


[ 14 ] 2. The pre-event image is subtracted from the 
sequence of CME event images. These calibrated images 
now display information on the excess (depletion) of light 
in units of mean solar brightness, MSB (i.e., a base- 
difference sequence). 

[ 15 ] 3. The excess (depletion) of MSB within each pixel 
in the image is converted into a number of electrons by 
using the Thomson scattering equations and by assuming 
all the electrons are located on a single plane determined 
by the fixed-phi J-map technique. 

[ 16 ] 4. The mass per pixel is calculated from each image 
by assuming a solar wind distribution of 90% H and 10% 
He. This corresponds to a mass of 1.97 x 10 -24 g per electron 
[Hildner et al, 1975]. 

[ 17 ] 5. The mass of the CME is then estimated by sum- 
ming up the values of all the pixels containing the CME. 

[is] Figure 2 is a combined and cropped image from 
the ST-A and includes a frame from COR-2, HI-1, and 
HI-2. Practically speaking, the pre-event image is usu- 
ally subtracted from the frames which occur later in time 
(i.e., when the CME has entered the field of view). How- 
ever, this frame may also be subtracted from images that 
occurred earlier. In Figure 2, the image from the HI-2A 
camera is one that is actually earlier in time than the 
pre-event image. Therefore, we are demonstrating the fea- 
sibility of using the pre-event image as a static background 
for frames earlier and later in time. 

[ 19 ] The pre-event image subtraction process is most 
suitable for relatively short time scales due to the steady 
state assumption. ITiis assumption therefore begins to 
breakdown for the long duration of the CME within the 
large field of view. Also, HI-2 is sensitive to the back- 
ground star field. This means that the ability to detect a 
CME motion is reduced if the star field is not appropri- 
ately removed during the initial image processing prior 
to making mass calculations. This paper is focusing on 
the potential to track the momentum flux and to monitor 
the relevant effects at Earth. For this reason, we choose to 
use the simplest image processing at this stage in order 
to emphasize the minimum capability for space weather 


forecasting. Further processing of the images should be 
able to improve the tracking of the CME momentum 
[Howard et al, 2012; Howard and DeForest, 2012] and may 
provide a more sensitive time series of mass flux at the 
LI point. 

[ 20 ] The majority of the current work on mass estimates 
uses data from coronagraphs which have a small field 
of view and assume the Thomson sphere is a flat plane. 
For the large field of view for both the HI cameras, this 
is inappropriate. When estimating the mass, the propa- 
gation direction, and therefore the angle away from the 
plane of sky (PoS), is important in calculating the amount 
of Thomson scattering from electrons, i.e., the pixel on the 
inner edge of the camera has a different angle away from 
the PoS to the outer edge. This modification to the calcu- 
lations is made in our work. In this paper, we treat each 
pixel individually as a different angle away from the plane 
of sky as measured along the line of sight. 

2.2. Semi-automated CME Mass 

[ 21 ] The calculation of the total CME mass is carried out 
by summing up the mass values of each pixel within the 
observed CME. This region can be defined in a variety 
of ways: 

[ 22 ] 1. The sector method [Vourlidas et al., 2010]. The 
observer manually defines a set of four boundaries that 
can be used for the entire sequence of images. These 
boundaries are defined between two position angles (PAs) 
and the inner and outer radial boundary. 

[ 23 ] 2. Region of interest (ROI) method [Vourlidas et al, 
2000; Subramanian and Vourlidas, 2007], The observer man- 
ually draws an outline of the CME. All the pixels within 
this ROI are considered to contribute toward the CME 
mass. This method calculates the total mass more accu- 
rately than the sector method but requires each frame 
to be considered individually and is therefore more time 
consuming. 

[ 24 ] 3. Graduated cylindrical shell model (GCS) method. 
This method uses the forward modeling technique 
developed by Themisien et al [2009] and Themisien [2011] 
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Figure 3. (top) The graduated cylindrical shell (GCS) 
model overlaid onto COR-2A and B mass images at 
15.24 UT, on 25 March 2011. (bottom) The measure- 
ments of the CME mass in COR-2 by using the GCS 
model as the CME boundary. 

to define the outer boundary of the CME. The different 
parameters in the model apart from the radial distance 
are manually chosen from comparing the model shape 
with both STEREO COR-2 images simultaneously. These 
parameters are then fixed and the radial distance is varied 
over the sequence of images as the CME propagates. The 
outer edge which is traced out by the model is used in the 
same manner as the ROI method defined above. This new 
GCS method is used in our paper. 

[ 25 ] The CME enters the COR-2A camera at 21.24 UT, 
on 24 March 2011, and approaches a mass of ~3 xlO 12 kilo- 
grams, which is of the order of magnitude of a typical CME 
[Vourlidas et al., 2010]. The mass estimates as measured 
from the GCS model are displayed in Figure 3 and show 
the typical increasing mass time profile due to the CME 
entering the field of view. This method enables the CME 
mass calculations to be made in a semi-automated fashion 
similar to the sector method, while using a more reliable 
ROI method for tracing the CME contour. In our example, 
the mass of the CME was slightly overestimated by <5% 
compared to using method 2. This is because we ensured 
that the entire CME structure was enclosed. For our event, 
which may not always be the case, it was relatively sim- 
ple as the CME had well-defined boundaries. However, 
as the CME propagated, the leading edge became slightly 
flattened in comparison to the idealized GCS model. This 
meant the nose of the GCS was progressively further into 
the heliosphere than the observations. This meant a few 
extra pixels of mass were included but less than would 
have been encountered by method 1. It was found that the 
accuracy of the mass measurements compared to method 


2 depended on ensuring that at least the entire area of 
the CME was included rather than minimizing the surplus 
area of background solar wind. 

[ 26 ] The GCS method used in this paper requires sev- 
eral parameters to be defined for each image. The results 
of propagation direction and flux rope axis were com- 
pared to the fixed-phi method (from remote observations) 
and constant-alpha force-free flux rope model (from in 
situ data; see section 3.1 for more details), respectively. 
The results were consistent with each other. The indepen- 
dently estimated propagation directions were within 10° 
of each other and the three independently estimated flux 
rope axes were within 30°. 

[ 27 ] Another option for measuring the mass flow is to 
investigate a fixed location in the heliosphere by defin- 
ing a small narrow rectangular box a few pixels wide (slit 
method). This method emulates a time series of mass over 
a fixed location such as a spacecraft at the LI position (see 
section 2.4). Currently, the minimal image processing used 
in this study does not allow for an accurate estimate to be 
made at LI. However, future studies using more advanced 
image processing and better tempo-spatial resolution as 
expected from the Solar Orbiter should provide the nec- 
essary data to advance the techniques in preparation for a 
possible mission to L5. 

2.3. Mass J-Maps 

[ 28 ] Originally developed for LASCO coronagraphic 
images, Sheeley et al. [1999] calculated that a small plasma 
packet moving at uniform speed would have an apparent 
acceleration and deceleration which is dependant on two 
variables: the radial velocity, Vr, and the angle of prop- 
agation between the CME-Sun-spacecraft, ft. Moreover, 
the observed acceleration profile (measured in time and 
elongation angle away from the Sun, a) is unique; there- 
fore by using an optimization routine [e.g., Savani et al., 
2009, 2010], an estimate of the propagation direction of 
the CME can be made. This technique (called fixed-phi 
method) has been shown to be much more useful over a 
large range of elongation angles [Williams et al., 2009; Davis 
et al, 2010] and produced very effective results in track- 
ing CMEs from the Sun to planetary systems where they 
were detected in situ [e.g., Davis et al., 2009; Ronillard et 
al., 2009; Mostl et al., 2009], As the premise on which this 
technique is built relies on the idea of a spatially narrow 
plasma packet and not a large three-dimensional (3-D) 
object traveling through the heliosphere, other attempts 
have been developed to mitigate against some of the sim- 
plified assumptions [Kahler and Webb, 2007; Lugaz et al., 
2009; Davies et al., 2012], 

[ 29 ] The propagation direction is estimated by mea- 
suring the elongation angle as a function of time. This 
is most conveniently done by tracking a single feature 
within a CME through the field of view of all the cam- 
eras by using a time-elongation map (J-map). As a plasma 
packet can be assumed to propagate radially away from 
the Sun, the radial cuts used in J-maps can be varied 
to suite the CME direction if it is away from the ecliptic 
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Figure 4. J-maps created for the position angle of the Earth (PA = 97°). (top) J-map dis- 
plays the results from total brightness images with the location of Earth and the manually 
selected track of the CME's leading edge, (bottom) A mass J-map created from a sequence of 
mass images. 


plane. In the case of our event, we choose a PA of 97° 
from ST-A, which corresponds to the direction to Earth. 
Until this paper, these J-maps have been created by pro- 
cessing MSB images into a running difference sequence, 
as shown in the top panel of Figure 4. The blue crosses 
are manually chosen from the image and are the points 
used in estimating the speed and propagation direction 
of the CME. Using the fixed-phi method, the radial speed 
and propagation direction were estimated as 371 kilome- 
ters per second and 80°, respectively. The bottom panel 
displays the same time-elongation tracks but created with 
a sequence of mass images. As can be seen between the 
two versions, the white-light J-map is clearer at identify- 
ing the propagation of the CME. This is partly because 
further smoothing and image processing were undertaken 
on the white-light J-map. It is hoped that further stud- 
ies into mass J-maps will be used to provide a time series 
estimate of the mass propagating over the Earth in a pro- 
cess similar to the slit method described above. This is 
because the slit method would effectively represent a hor- 
izontal line along a mass J-map. Although the tracks in 
the mass J-map displayed in our paper are "noisy," future 
studies may implement advanced image processing tech- 
niques that are currently under development [e.g., Howard 
et ah, 2012] to find interesting discoveries. 

[ 30 ] The tracking of the mass estimates can be used 
with a minimal number of assumptions about the expan- 
sion process to estimate the density of the CME. With 
improvements to the image processing, this informa- 
tion may be used as inputs to space weather forecasting 


models, instead of using LI in situ data. Although this will 
clearly be less reliable than LI data itself, it has the big 
advantage of being measured remotely and of the order of 
~48 hours in advance. This may prove to be a significant 
improvement for our forecasting capabilities. 

[ 31 ] Figure 5 displays the mass J-map along with the in 
situ measurements of density and velocity at LI. The track, 
which was made from the white-light J-map, is overplot- 
ted onto the mass J-map along with a dashed line to show 
the position of Earth. The track clearly intercepts the posi- 
tion of Earth at the same time as the CME-driven shock 
(and the associated density increase from the sheath) 
arrives at LI. The shock arrival at LI displays a sudden 
increase in momentum flux which then compresses the 
Earth's magnetosphere (see section 4). 

2.4. Mass Time Series 

[ 32 ] Under the premise that the CME momentum can 
be tracked and possibly be used as early solar wind input 
into space weather simulations, it is important to observe 
how a typical mass time series may look like and how 
it compares to the in situ density profile currently being 
used as simulation inputs. Figure 5 shows a normalized 
times series of the mass measurements (green dashed 
curve) taken at 3.8° elongation from ST-A, which is within 
the COR-2 field of view. The mass measurements were 
taken using the slit method at a plane of sky distance 
between 13.7 and 14 R$ for each frame. The data were then 
ballistically time shifted to the LI position by assuming 
the CME traveled at 370 kilometers per second and was 
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Figure 5. (top) The time series of the WIND proton density and velocity for the period 25 
March to 01 April 2011. A normalized mass time series from the mass images are overlaid 
onto the density axis (green). The mass time series is generated from the slit analysis at 3.8° 
elongation within the COR-2A camera. The time series is ballistically propagated to LI at 
a speed of 400 kilometers per second at the leading edge and 370 kilometers per second 
at the trailing edge, with a linearly decreasing speed profile. Qualitatively, the normalized 
mass time series shows a remarkable resemblance to the in situ density of the CME sheath, 
(bottom) The mass J-map with the manually selected CME leading edge. 


linearly expanding [Owens et al., 2005] so that the trailing 
edge of the mass measurements propagated 30 kilometers 
per second slower than the sheath leading edge (see also 
Figures 6 and 7). As our event initially propagates with a 
slow speed, we are able to assume the CME was swept into 
the solar wind and merely advected out to 1 AU [Siscoe and 
Schwenn, 2006]. However for faster CMEs, it is important 
to consider deceleration due to drag effects that change 
the arrival times at Earth [Gopalswamy et al, 2000; Cargill, 
2004]. The units of mass displayed have been normalized 
to suit the density profile. This has simply been carried 
out by dividing the measurements by 1.5 x 10 11 grams per 
cubic centimeter. A simple method of an expanding vol- 
ume during the CME's propagation is used to justify this 
value (see Appendix A for more details). 

[ 33 ] Currently only minimally processed images are 
being used and even with our simplified propaga- 
tion technique [e.g., Owens and Cargill, 2004; Lugaz and 
Kintner, 2012], and the qualitative profile of the mass mea- 
surements made from remote observations is remarkably 
similar to that observed in situ. It is clear that CMEs 
may undergo interactions during propagation to Earth by 
either solar wind distortions [e.g., Lugaz et al., 2008; Savani 
et al., 2010], deflections [Lugaz et al., 2012; Wood et al., 2012], 
or possible rotations [e.g., Shiota et al., 2005; Vourlidas et al., 
2011; Nieves-Chinchilla et al., 2012]. This could be the cause 
of the larger predicted mass measurements within the 


CME when compared to the in situ values. Or more sim- 
ply that the CME expanded more during the propagation 
and is more rarified by the time it reaches Earth. 

3. In Situ Observations 

[ 34 ] The shock associated with our case study event was 
detected at LI at 15.09 on 29 March 2011 with a jump in 
velocity of ~70 kilometers per second (the upstream speed 
was ~330 kilometers per second). Figure 6 shows the in 
situ parameters of the solar wind during the propagating 
interplanetary CME. The three vertical lines indicate the 
locations of the shock front (14.58 on 29), CME leading 
edge (23.39 on 29), and rear edge (09.52 on 31), respectively. 
The CME leading and rear edges were manually chosen by 
focusing on looking for a duration that includes a smooth 
rotation in magnetic field and a discontinuity on density. 
The focus on the smooth field rotation is in order to pro- 
duce reliable results from an in situ flux-rope fitting pro- 
cess. The magnetic field vectors are displayed in the RTN 
coordinate system such that the Normal (N) component 
is the measure of the magnetic field in the out-of-ecliptic 
direction (i.e., the north-south direction that is crucial in 
space weather predictions). This case study shows that 
the magnetic field in the CME and sheath region is pre- 
dominately northwardly directed and strongly so in the 
earliest half of the CME just behind the leading edge. This 
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Figure 6. In situ data from ACE and WIND at LI point. The first to sixth panels show the 
magnetic field in RTN Cartesian and spherical coordinates system, followed by proton veloc- 
ity, density, and temperature, respectively, in the seventh to ninth panels. The vertical lines 
from left to right display the positions of the CME shock and magnetic flux rope leading and 
trailing edge. 


CME topology was specially chosen for our analysis as it 
allows our study to isolate the space weather effects (e.g., 
potential strength of GICs) that are caused by the momen- 
tum flux and not from the resulting magnetic reconnection 
from a southwardly directed interplanetary magnetic field 
(IMF). 

3.1. Modeling Results for March 2011 

[ 35 ] Figure 6 displays the magnetic field profile of the 
optimized CAFF model (see Appendix B for details) within 
the first to sixth panels as black curves. The first to third 
panels represent the Cartesian vectors in the RTN coor- 
dinate system and the fourth to sixth panels represent 
the field vectors in the spherical coordinate system. The 
smooth rotation in the field is often more clearly seen in 
spherical coordinates while the importance of a southward 
Bz for predicting space weather events is better observed 
in the Cartesian coordinates. The orientation of the esti- 
mated flux rope axis direction is (0.4, -0.7, 0.6) in RTN and 


has a right-handed chirality. The optimized parameters for 
the axial magnetic field and the impact parameter are 14.0 
nanotesla and 0.1, respectively. This indicates the space- 
craft travel close to but slightly above the FR axis. The 
mean square error between the optimized model and the 
data, y, was <0.1, which represents a good fit to the data 
[Lynch et al., 2003]. The model fits were also compared to 
a non-force-free elliptical model [Hidalgo and Nieves-Ch- 
inchilla, 2012] to show consistent results (e.g., the axis 
orientation were within 30° of each other). 

3.2. Dst Index 

[36] There is a variety of ground responses at Earth 
from space weather disturbances, and measurements are 
often focused on different geographically localized pro- 
cesses/positions (e.g., due to local noon time, latitude, or 
ground resistivity). The disturbance time (Dst) index mea- 
sures the hourly values of the horizontal component of the 
Earth's magnetic field averaged from four near-equatorial 
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Figure 7. A comparison between the in situ data from 
LI (Bz component, velocity, and density) to the aver- 
aged global response by the Dst index. The normalized 
mass time series from the COR-2A is displayed on the 
density axis as in Figure 5. 


geomagnetic observatories. The Dst index has historically 
been used as an approximation of the global response 
to a space weather disturbance. The fluctuations in the 
Dst closely relate to the ring current and the other cur- 
rent systems (including the magnetotail current) within 
the terrestrial environment. The inverse proportionality 
relationship between the horizontal component of the 
magnetic field and the energy content of the ring current 
is known as the Dessler-Parker-Sckopke relation [Dessler 
and Parker, 1959; Sckopke, 1966]. 

[ 37 ] The stereotypical time profile of the Dst index dur- 
ing a geomagnetic storm driven by a fast CME displays an 
initial positive sharp rise (called the sudden storm com- 
mencement, SSC) which defines the arrival of the leading 
shock front onto the magnetosphere, a drop in value to 
zero (called the storm onset, SO), and then the main phase 
is characterized by a large negative decrease which rep- 
resents the period of strong southward Bz magnetic field. 
The sheath region between the leading edge of a CME 
and the shock front is often considered to be related to the 
period between the SSC and the SO (initial phase) during 
a geomagnetic storm. However, the phases of a geomag- 
netic storm can be time shifted under conditions where a 
strong southward Bz persists within the sheath region. 

[ 38 ] The case study CME analyzed in this paper was cho- 
sen to better understand the significance of its momentum 
flux as a driver of geomagnetic storms. As such, we chose 
to investigate a CME with a predominately northward Bz 
field. This allowed us to isolate the observed geomagnetic 
disturbance and assume the disturbance is solely due to 
the momentum flux and not due to reconnection between 
the CME and magnetosphere. For this reason, the usual 
main phase of a storm as seen in the Dst is not seen in 
our event (Figure 7). However, Figure 7 shows a signifi- 
cant rise in the Dst value and sharp fall during the initial 


phase, which represents the location of a sudden increase 
in momentum flux from the CME sheath region. The sud- 
den changes in the Dst show that the momentum flux 
of the CME leading edge is capable of making sudden 
changes to the ground magnetic field (i.e., cause a large 
dBldt). 

4. Terrestrial Response 

[ 39 ] Sophisticated MHD simulations in 3-D are becom- 
ing an increasingly effective tool for modeling solar wind 
transients such as CMEs and for predicting their geo- 
effectiveness at Earth. In this paper, we have employed the 
Space Weather Modeling Framework (SWMF) [Toth et al., 
2005] package which was executed at the Community 
Coordinated Modeling Center (CCMC) and operated at 
NASA Goddard Space Flight Center. The solar wind input 
data were chosen from the WIND spacecraft at LI and was 
ballistically mapped to the outer boundary of the BATS- 
R-US magnetospheric MHD model that was coupled to 
the Rice Convection Model in our simulations. The auroral 
conductances are driven by solar irradiance observations 
of F10.7 and field-aligned electric currents. 

[ 40 ] Figure 8 displays a 2-D cut of Earth's magneto- 
sphere (top) before the arrival of the CME-associated 
shock front and (bottom) when the sheath region of the 
CME is traveling over Earth's bow shock. The color table 
represents the density of the plasma and a few selected 
magnetic field lines are drawn to help distinguish between 
the terrestrial and heliospheric systems. The vectors show 
the solar wind direction. In the 3 hours between the frames 
shown in Figure 8, the bow shock is severely compressed 
from a location of approximately 16 Rs to —10 Rs ■ This 
compression significantly increased the density within the 
magnetosheath and was due to the sudden arrival of a 
larger momentum flux from within the sheath region of 
the CME. 

[ 41 ] Currently, this study uses the in situ data mea- 
sured at LI as the inputs into the SWMF simulations in 
order to predict realistic geomagnetic disturbances. How- 
ever, it is envisaged that with further development to the 
mass images, a reasonable density time series calculated 
remotely may be used as an input into the BATS-R-US 
code (see section 2.4). This may prove to be a valuable tool 
in improving any early warning systems by being able to 
provide an observationally predicted result that is at least 
—24 hours earlier than is currently possible. 

4.1. Observations of Magnetic Fluctuations 

[ 42 ] In order to better understand the geomagnetic 
effects driven by CMEs and in particular to predict the 
socioeconomic impacts from GICs, it is important to study 
the localized effects observed at specific locations on Earth 
during the arrival of the CME. In this paper, we study 
the momentum flux of CMEs and therefore choose to 
investigate the equatorially based ground stations. The 
higher latitude stations are likely to display magnetic field 
fluctuations that can be partly associated with auroral 
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Figure 8. Two-dimensional cut of the Earth's magneto- 
sphere and bow shock from the BATS-R-US simulation 
at CCMC. The color scale of the density is the same 
in both panels, (top) The shape and structure of the 
magnetosphere prior to the arrival of the CME shock, 
(bottom) The compressed magnetosphere at the time 
when the CME sheath is affecting the bow shock. 


magnetosphere-ionosphere dynamics, and are therefore 
not the focus of this paper. The shock arrival at LI (—15.00 
UT) defines the start of the geomagnetic disturbance. We 
investigate two locations on Earth that are determined to 
be locally noon and midnight at the shock arrival. The Vas- 
souras (VSS) station in Brazil is used as the locally noon 
station. The Kanoya (KNY) and Kakadu (KDU) stations in 
Japan and Australia, respectively, are studied as the locally 
midnight stations for our investigation. 

[ 43 ] In order to investigate the possible effects that a 
CME may have on GICs, the time derivative of the ground 
magnetic field (dB/df) must be studied. Figure 9 displays 
the dB/df at both local times and shows significant fluctu- 
ations at the arrival of the shock onto Earth's system. The 
vertical dashed line displays the time the shock arrived at 
LI and therefore the delay of —1 hour on the ground sta- 
tions predominately represents the propagation time of 
the solar wind between LI and the bow shock. The data 
are displayed with 1 minute temporal resolution in the 


Cartesian geographic (GEO) coordinate system provided 
by INTERMAGNET (www.intermagnet.org). 

4.2. Simulated dB/df 

[ 44 ] Quantifying and predicting ground magnetic field 
perturbations are vital to the space weather community. 
As such, the CCMC has developed a tool that is able to 
extract the ground magnetic field perturbations from the 
global MHD model outputs by integrating the results from 
the magnetospheric and ionospheric current systems. In 
particular, a summation of four separate current systems 
is used to make the predictions: (1) the current system 
in the magnetosphere (and magnetotail) above 2.5 R E , (2) 
field-aligned currents between 2.5 R E and 110 kilometers, 
(3) Hall current from the ionosphere, and (4) Pederson cur- 
rents in the ionosphere [Rastaetter et al., 2004; Pulkkinen 
et al, 2010b]. The performance and metric-based analy- 
ses between various modeling approaches were part of 
the Geospace Environment Modeling (GEM) 2008-2009 
challenge and is reported by Pulkkinen et al. [2011]. 

[ 45 ] Figure 9 shows the predicted estimates of the 
ground field as purple curves overplotted on the data. The 
simulated time series have been shifted earlier in time by 
15 minutes which we attribute to a small uncertainty in 
the input data which resulted from ballistically shifting 
the input LI data to the outer boundary of the simulation. 
We show that the simulated data also display a signifi- 
cant spike in the dB/df at the storm commencement which 
coincides with a similar magnitude to the "ground truth" 
observed by the magnetometer ground stations. Quali- 
tatively, the biggest limitation in the simulated results 
appears to occur in the inaccuracy of the negative dBzIdt 
component seen during the time derivative "spike." This 
could be due to the Bz component being especially sensi- 
tive to geomagnetic induction effects which are not taken 
into account in our simulation. 

[ 46 ] It is worth noting that only middle- and high- 
latitude magnetometer stations were included in the ear- 
lier GEM challenges even though the Dsf index addresses 
the low-latitude disturbances. This is due to the global 
MHD approach being implemented in the first-principle 
models; they could only capture the ionospheric output 
at high latitudes by using the Biot-Savart law to inte- 
grate over the ionospheric electric currents system. In 
this paper, we have alleviated such a constraint by cou- 
pling the global MHD models to the inner magnetospheric 
models and thereby capturing the ring current dynamics 
and magnetospheric current systems, therefore provid- 
ing the required ionospheric response at low latitudes 
[Yu et al., 2010]. 

4.3. Geomagnetically Induced Currents, GICs 

[ 47 ] Reliable estimates of GICs and the geoelectric field 
require accurate knowledge of the local geological con- 
ditions as well as the dB/df. As the global distribution 
of the conductivity from the surface to the upper man- 
tle (depths of several hundred kilometers) is not well 
known, estimating GICs cannot be arbitrarily made for 
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Figure 9. The time derivative of the magnetic field measured at the low-latitude Vas- 
souras (Brazil) and Kakadu (Australia) ground stations. These stations represent the local 
noon and midnight of the CME shock front onto the bow shock, respectively. The purple 
curves overlaid onto the data are the estimated values at these locations from the CCMC 
simulations. 


any location on Earth. However, for local environments 
that have historically been susceptible to GIC events [e.g., 
Pulkkinen et ah, 2005; Ngwira et al., 2008; Pulkkinen et ah, 
2010a; Torta et al., 2012], the ground structure is known 
and therefore this study could be replicated for larger geo- 
effective CME events and early warning predictions can 
be made as to their level of susceptibility. It is envisaged 
that as the ground structure for more locations around 
the globe become recognized, the framework presented in 
this paper will provide the necessary steps to improve GIC 
forecasting. Of course, a simple approach may be created 
to estimate the extrema of possible GICs by using real- 
istic extreme ends of the conducting (British Columbia, 
Canada) and resistive (Quebec, Canada) ground struc- 
tures [Pulkkinen et al, 2008]. But detailed investigation into 
these possibilities goes beyond the scope of this work. 


5. Discussion 

[ 48 ] Currently, observations of the solar magnetic fields 

are used as inputs for the background solar wind when 

simulating the heliosphere for space weather predictions. 
However, the CME itself is usually set to an approximate 

and generic size. Details of the CME that is included in the 
magnetospheric simulations are only provided by mea- 

surements made in situ at LI, which is ~1 hour before its 

arrival at the Earth's bow shock. By using remote observa- 
tions from coronagraphs (~15 Rs), observational estimates 
of the CME can be made for both the Enlil heliospheric 
simulation and as early initial-attempt inputs to BATS-R- 
US magnetospheric simulations. The use of remote obser- 
vations means that BATS-R-US simulations can be carried 
out ~48 hours prior to waiting for the CME to propagate 
to LI. 


[ 49 ] The focus of this paper was to estimate the CME 
mass as a time series close to the Sun, which was then 
ballistically propagated to LI. While previous studies 
have been carried out to better estimate the arrival times 
of a CME [Owens and Cargill, 2004; Taktakishvili et al., 
2009, 2011], extensive plane-of-sky speed measurements 
have been made over a large number of CMEs [St. Cyr 
et al., 2000; Yashiro et al., 2004]. These studies could be 
used to remotely estimate a velocity time series which 
can then be used to propagate the CME to LI. The 
ability to remotely estimate both the velocity and mass 
allows early prediction of the momentum flux arriving to 
Earth. 

[ 50 ] Figures 5 and 7 display the remotely observed mass 
time series which have been artificially normalized to 
suit the in situ number density. The mass time series 
were arbitrarily divided by 1.5 x 10 5 * * * * * 11 g in order to dis- 
play a qualitative profile of the same order of magni- 
tude. Clearly, if these mass time series are to be used 
as inputs to early space weather simulations, then an 
appropriate method to scale the mass estimates to a den- 
sity value is required. Three basic methods can be used 
as follows: (1) defining the CME's 3-D volume by, for 
example, using the GCS model [Themisien et al., 2006], 
this volume can then be radially propagated out to 1 
AU where the density can be estimated; (2) using the 
empirical formulas proposed by Vourlidas et al. [2010]; 
and (3) using a combination of white-light images and 
off-limb spectroscopy to directly measure the densities 
across the CME body. Future investigations may attempt 
to solve the most appropriate method with more events. 
Method 1 has been used in Appendix A to justify the 
order of magnitude used in this paper. With appropri- 
ate density and velocity estimates made remotely to first 
approximation, the results can be adjusted in an ensemble 
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method for making space weather predictions along 
with variety of magnetic field estimates [e.g., Lin et al., 
2000; Gopalswamy and Yashiro, 2011; Savani et al., 2012b]. 

6. Conclusions 

[ 51 ] In this paper, we show the first results for produc- 
ing density measurements from remote observations that 
are comparable to their equivalent in situ time series. 
Therefore, this paper for the first time shows that a time 
series of data can be estimated remotely and be used to 
make forecasts of GICs at Earth, while a previous attempt 
by Pulkkinen et al. [2009] to forecast GICs from remote 
observations only used a generic pressure pulse for a 
CME within the Enlil model. We also confirm that the 
compression of the magnetosphere from the momentum 
flux of a CME is a significant variable in predicting geo- 
magnetically induced currents. Our results are capable of 
producing qualitatively reliable estimates of the densities 
upstream of Earth's bow shock. These results indicate the 
possibility of using remote observations at a heliocentric 
distance of ~15 R s to estimate the solar wind density pro- 
file arriving at Earth and show that these estimates can 
be used as part of a preliminary early warning system for 
space weather predictions. 

[ 52 ] This paper focuses on a case study event of a CME 
which displays a strong northward Bz component in the 
magnetic field rather than a fast geo-effective event. The 
CME propagated in the inner heliosphere between 25 and 
31 March 2011 and is observed in both ST-A and ST-B. The 
event allowed our study to attribute all the ground-based 
effects at low latitudes to the compression of the magneto- 
sphere and not to the magnetic reconnection in the case of 
a southward directed Bz. The compression of the magne- 
tosphere is found to be the result of the larger momentum 
flux (larger density and velocity) in the sheath region of 
the CME. 

[ 53 ] Previous studies have estimated the mass of a CME 
as it traveled through the field of view either by man- 
ually measuring the region of interest (ROI) around the 
CME for each individual frame or by estimating it within 
a fixed sector and thereby assuming that the inclusion of 
the additional solar wind mass is not significant in relation 
to the CME mass. Here, we improve the process by semi 
automating the ROI while refraining from the manual 
selections within each frame. We do this by using the grad- 
uated cylindrical shell model [ Themisien , 2011] to define 
the ROI and manually fixing the parameters by eye at the 
beginning and then by only varying the radial distance for 
each consecutive frame. 

[ 54 ] To track the CME mass from the remote obser- 
vations to LI, we employed the J-map technique for 
a sequence of mass images in the STEREO data. The 
J-map technique has previously been extensively used 
with white-light images to estimate the propagation direc- 
tion and arrival time of a CME [e.g., Sheeley et al., 1999; 
Rouillard et al., 2009]. However, we show that the mass 
J-map technique can be used to highlight the mass 


intensities traveling toward Earth. As mass images are 
estimated from base-difference frames rather than run- 
ning differences, the background star field becomes sig- 
nificant in the HI-2 field of view. In this paper, we have 
chosen to use a minimal amount of image processing in 
order to estimate a base confidence level for the uncer- 
tainties in the mass values. Future studies may consider 
a comparison between the minimal processed images to 
more advanced techniques. 

[ 55 ] As the mass images in the vicinity of Earth (HI-2 
field of view) was not clearly visible in the images due 
to the interference of the background star field, we chose 
to estimate the mass as a time series at a fixed heliocen- 
tric location within the COR-2A field of view (plane of sky 
distance of 14 Rs). The mass time series was then time 
shifted by assuming a leading edge speed of 400 kilome- 
ters per second and having a linear expansion profile such 
that the trailing edge of the mass profile was traveling at 
370 kilometers per second. The qualitative comparison of 
our estimated mass profile to the in situ density measure- 
ments at LI is remarkably similar for the sheath region 
between the shock front and leading edge. Future studies 
to investigate CMEs may not necessarily have the capabil- 
ities to track the CME along the entire Sun-Earth line. As 
such, our ballistic propagation approach for the mass mea- 
surements from ~15 Rs to LI may serve as a useful tool for 
estimating a time series for the mass. 


Appendix A: Density Normalization 

[ 56 ] In this paper, we have chosen to artificially divide 
the time series of the mass measurements by a constant 
value of 1.5 x 10 n g. This was used to convert a mass 
times series that was estimated remotely into a density 
time series which was later compared with the in situ mea- 
surements. Here we carry out some preliminary work to 
justify the number used with a simple volume expansion 
method. The work below is intended to provide an order 
of magnitude justification; however, further work would 
benefit from a more detailed approach as suggested in 
section 5. 

[ 57 ] We assume that a small volume (defined by the 
size of a few pixels in the COR-2A camera) measured 
at ~10 Rs has a shape of a cylinder, such that the cir- 
cular cross section is within the plane of sky and the 
length (W 0 ) is defined along the line of sight. First, we 
assume that the length of the cylinder can be estimated 
from the GCS model (see Figure Al for details). At 15.24 
UT on 25 March 2011, the height of the legs ( h ) and the 
half-angular width (a) from GCS as defined by Thernisien 
et al. [2006] were estimated as 9.78 Rs and 32.5°, respec- 
tively. With trigonometry, we calculate that the cylindrical 
length should be 10.51 Rs- For this time, we noted that 
the plane of sky position of the CME nose was 10.5 Rs- 
The same process was carried out for an image at 17.24 
UT to deduce that the heliocentric distance of the CME 
was to first approximation equal to the cylindrical length. 
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Figure Al. (a) The parameters from Tliemisien et al. [2006] that are used to estimate the 
length of the CME along the line of sight, (b) The initial volume element defined to be 
observed with the COR-2A camera, (c) The final dimensions of the volume element at LI, 
where the in situ measurements are taken. 


Wq. Therefore, at LI where the in situ measurements are 
made, we estimate the cylindrical length as W\ = 210 Rg. 

[ 58 ] To estimate the cross-sectional area, we assume the 
cross section expands at a uniform speed of V ex = 30 kilo- 
meters per second (which we estimated from the in situ 
measurements). The total time of expansion, t, is the same 
as the propagation time and dependent on the bulk flow 
speed which we assume to be 370 kilometers per second. 
Assuming the propagation distance, R = 200 Rs, we find 
that the radius of the circular cross section at LI is 16.22 R$ 
and follows 



[ 59 ] Therefore, our final volume element (I/) can be 
estimated by 

i >-' i (£) 2r2,v '' (A2) 

to give 5.85 x 10 37 cubic centimeters. 

[ 60 ] It therefore follows that a first approximation for a 
normalization constant (at 0 ) is 


[61] From the constants used in our example, k 0 = 1.02 x 
1 0 1 1 y -1 cm -3 . This is similar to the normalization value of 
1.5 x 10 11 used in this paper and is certainly of the cor- 
rect order of magnitude. This small discrepancy might be 
solved with more advanced calculations. 

Appendix B: Flux Rope Fitting 

[ 62 ] Simple determination of parameters from a model 
fitting procedure is one of the best ways to quickly esti- 
mate the global properties of an interplanetary CME 
(ICME) from in situ measurements. The first model 
to be successfully optimized to magnetic clouds (MCs) 
was a constant a, force-free (CAFF) flux rope model by 


Lepping et al [1990]. Since then, several other attempts 
have been made to improve the results between observa- 
tions and models [Owens et al, 2006; Mulligan and Russell, 
2001; Hidalgo et al, 2002; Hidalgo and Nieves-Chinchilla, 
2012; Owens et al, 2012]. We use a simplified modifica- 
tion of the Lepping et al. [1990] model. The magnetic field 
of the model is described by the Cartesian form of the 
Bessel function for a force-free flux rope. In fitting the 
flux rope model to the data, the first step is to estimate 
the orientation of the cylindrical axis. Then by determin- 
ing the chirality and varying both the distance of closest 
approach and axial magnetic field strength as free param- 
eters within a computational code, we are able to make 
predictions of ICME parameters. 

[ 63 ] The method employed to determine the MC axis 
direction involves a technique called minimum variance 
analysis (MVA) [Sonnerup and Cahill, 1967; Steed et al, 
2008]. MVA was originally developed to determine the 
local normal for a tangential discontinuity (TD) in the 
magnetic field. This method calculates the orthogonal set 
of vectors in which the variance of the magnetic field in 
question is at minimum (el), maximum (e3), and inter- 
mediate (e2). For each of these eigenvectors, there is a 
corresponding eigenvalue, often represented as A 3 , A 3 , and 
A 2 , respectively. 

[ 64 ] For the purposes of investigating TDs, the mini- 
mum eigenvector calculated provides the normal to the 
plane. Flere the ratios of eigenvalues (Ai) and the other two 
are of concern as they indicate the reliability in the nor- 
mal. For investigating the axial direction of a MC, we must 
concern ourselves with the intermediate direction, with 
MVA applied to the time period that contains the flux rope 
only. The axis of the MC therefore lies in the intermediate 
direction, but the distinction of the field variance along the 
coordinate axes reduce as the spacecraft trajectory moves 
further from the central MC axis. 

[ 65 ] When investigating the intermediate direction, the 
eigenvalue ratios between A 2 /Ai and A 3 /A 2 are of con- 
cern. Ideally, both ratios should be > 10, but recent work 
using MVA has shown that ratios greater than 2 are ade- 
quate in ICME context [Bothmer and Schwenn, 1998]. If 
the ratio of either A 2 /A 3 or A 3 /A 2 is small, then the axial 
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direction becomes ambiguous within the plane of the two 
eigenvectors. 

[66] The scenario above is described for the ideal case, 
where the spacecraft passes through the center of the flux 
rope. However, this is usually not true. The distance of 
closest approach, Y 0 , is an important parameter in deter- 
mining the predicted field profile. It is defined in dimen- 
sionless units as the closest approach distance divided by 
the radius of the flux rope. Thus, Y 0 varies between 1 and 
-1 for a spacecraft traveling above and below the MC cen- 
ter, respectively. Rees [2003] simulated idealized flux ropes 
with spacecraft trajectories away from the center and then 
carried out MVA over the time profiles. He found that 
the error in the axis angle increased with an increasing 
|Y 0 |. The results suggest that an impact parameter of >0.5 
should be dealt with caution as the error may be of the 
order of ~15°. The uncertainty in ascertaining the axis 
orientation [Al-Haddad et al, 2013] is one of the reasons 
why more complex models choose to begin their analy- 
sis with MVA and then later find the optimal direction by 
introducing two extra free variables. 

[67] High-frequency noise is a source of error for MVA 
in the interplanetary CME context. MVA analyzes the vari- 
ance in the magnetic field, and therefore the large noise 
fluctuations can yield incorrect eigenvectors. To limit this 
effect, 10 or 15 minute averaged data are often used when 
fitting a modeled rope to the data. It has been shown 
that for good events, the time resolution between 10 and 
60 minutes has a negligible effect on the axis orienta- 
tion [Lepping et al, 2003]. Subtly, the lower resolutions are 
often used for large events, with a reason to maintain 
the total number of data points within the MC to about 
N = 40 [Lepping et al, 2006]. Although this is not essential, 
it is desirable when comparing various cases against each 
other in a statistical manner. 

[68] Once the axis of the model MC is determined, the 
in situ magnetic field data are rotated into the frame of the 
MVA axes. A sequence of flux ropes is then created and 
compared to the data by analyzing the mean square error 
(/) coefficient. The / coefficient measures how well each 
rope fits to the data and is defined in the same manner as 
Lynch et al [2003]. 

[ 69 ] The time interval between the start and end of the 
MC is identified manually by eye. This defines the data 
range which is compared to the model and fixes the size of 
the MC model. By using the distance of closest approach, 
the axial field strength, and chirality (a concept of rotation 
in the magnetic field, either clockwise or anticlockwise) as 
free parameters, various simulated MCs can be generated. 
Each simulated MC is compared to the data by measur- 
ing the x coefficient. The simulated MC with the smallest 
X value is regarded as the best fit result. The free vari- 
ables that create the model are varied by using a downhill 
simplex method developed by Nelder and Mead [1965]; this 
is a nonlinear optimization routine designed to minimize 
the x coefficients. This approach is faster than a stan- 
dard grid search and can be easily manipulated to deter- 
mine the free variables to larger significant figure. The 


optimized free parameters are the outputs of the model. 
These results are used again to re-create the optimized 
model. The field vectors from this model are then rotated 
into the spacecraft frame and plotted on top of the 
observed data for visual confirmation. 

[ 70 ] It is worth noting here that the CAFF model 
described above assumes a static ICME. That is to say, the 
model plotted shows a time series obtained by taking a 
radial cut through a flux rope at a fixed time. Many ICMEs 
observed at terrestrial distances take —24 hours to tran- 
sit over a spacecraft; therefore, it is important to note that 
ICME seen at the end of the transit has evolved and is dif- 
ferent from the same object at the beginning of the transit. 
A more accurate representation would be to include a time 
series of a flux rope evolving in time past a fixed point in 
space. 
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